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ABSTRACT 

Flow  diagrams  are  presented  for  a  computer  program  which  calculates  the  coordinates 
in  space  of  atoms  of  a  chain  molecule,  for  all  possible  rotational  isomers  consistent  with 
given  bond  lengths, bond  angles,  and  restrictions  on  rotation  about  single  bonds.  A  Fortran 
program  illustrates  one  form  of  this  computation,  and  modifications  are  outlined  for  appli- 
cation to  other  than  chain  molecules. 
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FLOW  DIAGRAMS  FOR  COMPUTATION  OF  SPATIAL  LOCATIONS  OF  ATOMS 
OF  A  CHAIN  MOLECULE  UNDERGOING  INTRAMOLECULAR  ROTATIONS 

C.    Roland  Eddy 


This  bulletin  supplements  the  note  that  appeared  in  the  Journal  of  Chemical  Physics  (4)* 
by  giving  flow  diagrams  for  a  digital-computer  program. 

The  equations  on  the  flow  diagrams  are  obtained  from  the  equations  in  the  Journal  by 
two  transformations:  (1)  The  subscript  i  is  given  in  turn  the  three  values  0,1,  2  to  cover 
the  three  coordinates.  Arithmetical  operations  on  i  are  performed  modulo  3.  (2)  In  order 
to  make  the  computation  iterative,  an  integral  number  p  is  assigned  to  each  atom  in  the 
chain,  beginning  at  one  end,  where  the  first  three  atoms  are  placed  in  known  positions  to 
specify  the  coordinate  system.  For  any  value  of  p  during  the  course  of  the  iteration,  Atom 
B  is  designated  Atom  p,  Atom  A  is  p-1,  Atom  C  is  p+1,  etc.  Direction  cosines  Sp,  tp,  Up, 
are  for  the  directed  line  from  Atom  p-1  pointing  toward  Atom  p.  Bond  length  Lp  is  the  dis- 
tance from  Atom  p-1  to  Atom  p.  Bond  angle  Vp  is  for  the  angle  at  Atom  p,  expressed  as  the 
supplement  of  the  angle  formed  by  the  three  atoms  p-1,  p,  p+1.  Dihedral  angle  Jp  is  the 
angle  of  rotation  of  the  BCD  plane  about  the  bond  joining  Atoms  p  and  p+1.  The  zero  of  J„ 
is  taken  as  the  trans  position  of  Atom  p+2  relative  to  Atom  p-1,  with  all  four  atoms  p-1,  p, 
p+1,  p+2  lying  in  the  same  plane.  Positive  Jp  is  defined  as  counterclockwise  rotation  of 
Atom  p+2  from  the  Jp=0  position,  as  viewed  from  Atom  p  looking  toward  Atom  p+1.  These 
designations  of  the  variables  are  illustrated  in  figure  1. 


Fig.    1 
Directed  lines     and  angles 
for  any  given  portion  of  the 
molecular  chain  as  a  func  - 
tion  of  subscript  p. 


Figures  in  parentheses  refer  to  Literature  Cited  at  end  of  report. 


For  starting  the  iterative  procedure,  it  is  necessary  to  specify  numerical  values  for  the 
coordinates  of  Atoms  0  and  1,  along  with  their  direction  cosines.  Direction  cosines  point- 
ing toward  Atom  0  (s0,  t0,  uQ)  have  no  molecular  significance.  They  merely  control  the 
plane  in  which  Atom  2  is  to  be  placed,  and,  in  so  doing  they  define  the  orientation  of  the 
coordinate  frame  of  reference.  The  program  automatically  sets  Vo=90°and  Jo=0°;  therefore 
s_,  tQ,  u0  must  be  chosen  as  a  line  perpendicular  to  the  bond  between  Atoms  0  and  1,  point- 
ing toward  Atom  0  from  a  position  trans  to  Atom  2,  and  lying  in  the  plane  desired  for  Atoms 
0,  1,  and  2.  Interrelations  within  the  plane  chosen  for  Atoms  0,  1,  and  2  are  illustrated  in 
figure  2. 

ATOM     2 


Fig.  2 
Directed  lines  in  the  plane  of 
Atoms  0,  1,  and  2,  showing 
the  relation  between  sQ,  tQ,  u0 
and  Atom  2  for  starting  the  it- 
erations. 
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As  p  is  incremented  along  the  chain,  the  coordinates  of  each  atom  are  computed  by  the 
basic  computational  unit  and  printed  out,  until  the  other  end  of  the  chain  has  been  reached  . 
This  gives  the  conformation  of  the  molecule  for  one  starting  set  of  values  of  the  J's; for  ex- 
ample,   those  corresponding  to  a  fully  extended  chain. 


The  index  p  then  follows  either  of  two  different  sequences,  depending  on  whether  it  is 
incrementing  or  decrementing.  During  the  incrementing  sequence,  for  each  new  value  of 
p,  Jp  is  reset  to  its  lowest  value.  During  the  decrementing  sequence,  for  each  new  value  of 
p,  the  existing  value  of  Jp  is  increased  to  its  next  higher  value  and  tested.  If  Jp  has  already 
passed  through  all  of  its  allowed  values,  p  continues  to  decrement;  if  Jp  has  not  yet  passed 
through  all  of  its  allowed  values,  p  is  switched  to  its  incrementing  sequence.  Whenever  a 
Jp  is  changed  to  a  new  allowable  value,  the  location  of  the  corresponding  D  atom  (numbered 
p+2)  is  computed.  Whenever  p  exceeds  its  upper  limit,  it  is  switched  to  the  decrementing 
sequence.  When  p  falls  below  its  lower  limit,  computation  ceases  because  all  J's  have 
passed  through  all  their  allowed  values  and  the  atom  locations  for  all  possible  rotational 
isomers  have  been  printed  out. 


The  integer  n  specifies  the  upper  limit  of  p.  Since  the  chain  of  atoms  runs  from  Atom  0 
to  Atom  n+2,  n  must  be  assigned  a  value  lower  by  3  than  the  number  of  atoms  in  the  chain  . 
For  iterating  J's,  the  initial  angles,  Jpo,  must  be  specified,  along  with  the  incremental 
angle  for  each  bond,  AJp,  based  on  rotational  restrictions  of  the  assumed  model.  The  up  - 
per  limit  of  Jp  is  taken  as  360°.  If  computation  for  only  one  rotational  isomer  is  desired, 
setting  each  AJp  equal  to  361°  will  cause  the  computer  to  stop  after  passing  along  the  chain 
only  once. 


Cases  that  are  mathematically  possible  but  physically  impossible  are  eliminated  by  as  - 
signing  a  minimum  radius  to  each  atom,  Rp.  As  each  Datom  is  located,  its  distances  from 
all  previous  atoms  are  computed.  If  any  one  of  these  distances  is  less  than  the  sum  of  the 
two  corresponding  minimum  radii,  that  rotational  isomer  is  eliminated  from  further  com  - 
putation.   If  the  chain  atom  carries  no  attached  atoms,  its  van  der  Waals  radius  can  be  used 


for  Rp.     If  the   chain  atom  carries   attached  atoms,   the  effective  radius  of  the  whole  group 
can  be  used.    For  example,    2.  0  Angstroms  can  be  used  for  methyl  or  methylene  groups  (9). 

Figures  3  and  4  give  flow  diagrams  for  the  case  of  a  series  of  atoms  connected  as  a 
single  chain.  This  computation  has  been  programmed  for  Bendix  G-15D  and  IBM- 1620  com- 
puters. *  The  program  is  illustrated  in  more  detail  in  Fortran  II  formalism  (pages  9  -  11). 
Included  is  an  additional  feature  omitted  from  the  flow  diagram  for  clarity:  The  Fortran 
program  places  up  to  two  side  atoms  on  each  chain  atom,  such  as  hydrogens,  halogens,  car- 
bonyl  oxygen,  etc.  This  program  is  not  necessarily  the  most  efficient  way  in  which  this 
computation  can  be  expressed  in  Fortran  form  but  is  chosen  to  illustrate  the  principles  used 

Specifications  for  this  program  are  given  on  page  12.  Figure  5  (page  8)showsthe  geometry 
of  the  main  chain  in  the  symbolism  of  the  Fortran  program.  Figure  6  shows  the  geometry 
for  the  side  atoms.  Table  1  (page  14)  gives  correlations  among  the  different  systems  of 
symbols,    made  necessary  by  the  Fortran  limitations  on  integers  and  subscripts. 

EXAMPLE 

As  an  example  of  the  use  of  this  program,  it  was  applied  to  homoallylic  alcohols  of  the 
general  type 


r 


R  —  CH=CH-  CH~2—  CH—  O—  H 
12  3  4 

to  determine  how  many  rotational  isomers  bring  the  OH  near  enough  to  the  double  bond  to 
permit  intramolecular  hydrogen  bonding.  For  standard  bond  lengths  and  angles,  the  fol- 
lowing values  were  used  (1,  2,  3,  6,  7,  10):  0^=02  1.337,  C2-C3  1.506,  C3-C4  1.533,  C4-O 
1.43,  O-H  0.960,^_C-C-C  112.5°,  /_OC-C  123°,  /_C-C-0  tetrahedral,  /_C-0-H  108°.  For  J 
angles  of  rotation  about  bonds  C3-C4  and  C4-O,  values  of  0°,  120°,  240°,  were  used,  corres- 
ponding to  the  accepted  potential  minima.  Since  the  potential  barriers  around  a  single  bond 
adjacent  to  a  double  bond  are  lower  and  less  definitely  located  (5,  8)  than  for  purely  single 
bond  systems,    12  equidistant  values  of  J  about  bond  C2-C3  were  assumed  beginning  with  0°. 

Of  the  resulting  108  rotational  isomers,  the  computation  finds  18  isomers  having  the 
hydroxyl  hydrogen  within  3.  0  A  of  the  center  of  the  double  bond  (chosen  as  the  origin  of  co- 
ordinates) and  10  isomers  with  the  hydrogen  less  than  2.5  A  from  center  of  double  bond. 
Table  2,  (page  15)  gives  coordinates  for  the  12  isomers  whose  hydroxyl  groups  are  most  fa- 
vorably placed  for  intramolecular  hydrogen  bonding  when  standard  values  of  angles  and  bond 
lengths  are  used.  Influence  of  nonstandard  angles  on  these  coordinates,  with  correlation  to 
infrared  spectra,   will  be  discussed  in  another  publication. 

MODIFICATIONS 

Other  programs   could  be  written  around  the  basic  computational  unit  to  accommodate 
modifications  such  as  the  following: 

(1)    Instead  of  assigning  the   coordinate   system   at  the  beginning  of  the  chain,  the   co- 
ordinate system  could  be  defined  by  assigning  locations  of  any  3  consecutive  atoms     anywhere 

*  Mention  of  a  proprietary  product  in  this  publication  does  not  constitute  a  guarantee  or 
warranty  of  the  product  by  the  U.  S.  Department  of  Agriculture  and  does  not  imply  its  ap- 
proval by  the  Department  to  the  exclusion  of  other  products  that  may  also  be  suitable. 


Figure  3 
Flow  diagram,    Part  1,    showing  basic  computational  unit  and  initial  steps. 
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\Print   out:    p+2,    x^ 2,yp+2  >  2p+2>    Jp     / 


Figure  4 
Flow  diagram,  Part  2,  showing  recycling  sequence  for  case  of  skeletal  atoms  of  chain  molecule 
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Decrementing  Sequence 
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Return  to  Basic  Computational  Unit,  Beginning  at  Block  2. 


in  the  molecule.     Then  chains  would  be  computed  in  both  directions  from  this  point. 

(2)  Branching  of  chains  could  be  accommodated  by  minor  changes  in  the  recycling  pro- 
cedure. 

(3)  To  accommodate  a  ring,  chains  could  be  computed  in  both  directions  around  the 
ring  until  they  meet  in  one  atom.  The  program  would  then  select  only  those  rotational  iso- 
mers which  give  two  sets  of  identical  coordinates  for  this  one  atom.  If  necessary ,  the 
computer  could  be  programmed  to  adjust  one  or  more  of  the  Jangles  by  trial  and  error  until 
the  two  locations  of  the  common  atom  become  identical.  Similarly  the  V's  and  L's  could  be 
adjusted  if  required.  Trial-and-error  methods  are  not  necessary  when  parameter  values 
are  previously  known  that  permit  a  closed  ring. 

(4)  For  each  rotational  isomer,  any  desired  property  of  atomic  locations  or  of  direc- 
tion cosines  could  be  computed,  such  as  electric  moment,  interaction  energy,  or  moments 
of  inertia.  By  making  assumptions  of  the  population  of  each  rotational  isomer  in  a  sample 
of  the  compound,  average  values  of  these  properties  could  be  computed  for  correlation  with 
experimental  measurements. 

(5)  One  or  more  of  the  angles  or  bond  lengths  could  be  varied  in  order  to  determine 
their  effects  on  molecular  properties  or  on  the  number  of  rotational  isomers  permitted  by 
the  requirements  on  minimum  interatomic  distances. 

(6)  The  program  can  be  very  conveniently  used  for  generating  data  for  the  orthogonal 
projections  of  Vandenheuvel  (11).  It  is  only  necessary  to  plot  two  of  the  coordinates  of  each 
atom  and  ignore  the  third  coordinate. 


Figure  5 

Geometry  of  Fortran  program 
for  main  chain 


BF  is  perpendicular  to  both  AB  and  BG 
BG  and  CG'  are  parallel  to  each     other 

and  both  are  perpendicular  to  BC. 
H2  lies  in  the  plane  of  B,  C,    and  D     and 

is  perpendicular  to  both  CG'  and  BC. 


Figure  6 


Geometry  of  Fortran  program  for  side  atoms 


CG'  is  perpendicular  to  the  plane  of  B,C,  and  D. 
H?  lies  in  the  plane  of  B,    C,    and  D  and  is  per- 
pendicular to  BC. 
Ho  is  perpendicular  to  CG'. 
-H3  bisects  the  angle  BCD. 

H4  lies  in  the  plane  of  H- and  CG', making  angle 
6  with  H3. 


Fortran  II  Program  for  Computing  Backbone  Chain 
Plus  Two  Side  Atoms  on  Each  Chain  Atom 

*1005 

DIMENSION  A(3,22)  ,B(3)  ,E(22)  ,F(3)  ,G(3)  ,UJ(20)  ,UO(20)  ,DU(20)  , 
1R(21)  ,V(21)  ,X(3,  22)  ,H2(3)  ,H3(3)  ,H4(3)  ,EH(21)  ,HCH(21)  ,XLINE(3)  , 
1XH(3) 
RADIN=1.745329252E-02 
101  READ  200,N,NSMPL 

200  FORMAT  (12,14) 
DO  1  L=l,2 

READ  201,X(1,L),X(2,L),X(3,L) 

1  READ  201,A(1,L),A(2,L),A(3,L) 

201  FORMAT  (3F14.10) 
KK=N+3 

DO  2  L=1,KK 

READ  202 , E (L) , V (L) , UO (L) , DU (L) , R (L) 

READ  202 , HCH (L) , EH (L) 

202  FORMAT  (5F14.10) 
V(L)=V(L)*RADIN 
UO(L)=UO(L)*RADIN 
DU(L)=DU(L)*RADIN 

2  HCH(L)=HCH(L)*RADIN 
L=l 

LF=1 
LD=1 

UJ(1)=0.0 
VSIN=1.0 
ISOMR=l 
GO  TO  503 
500  LF=L 

51  IF(V(LF))53,52,53 

52  LF=LF-1 
IF(LF-1)54,54,51 

54  VSIN=1.0 
GO  TO  503 

53  VSIN=SINF(V(LF)) 
503  LX=L+1 

DO  3  1=1,3 

3  B(I)=A(I,L+1) 
IF(VSIN)50,504,50 

50  DO  4  1=1,3 

GO  TO  (5,6, 7), I 

5  J=2 
K=3 

GO  TO  4 

6  J=3 
K=l 

GO  TO  4 

7  J=l 
K=2 

4  F(I)  =  (A(J,LF)*A(K,LF+1)-A(J,LF+1)*A(K,LF))/VSIN 


Fortran  II  Program  (Cont.  ) 

504  LY=L+2 

505  DO  8  1=1,3 
GO  TO  (9, 10, 11), I 

9  J=2 
K=3 
GO  TO  8 

10  J=3 
K=l 
GO  TO  8 

11  J=l 
K=2 

8  G (I) =F (I)*COSF (UJ (L) )+ (F (J)*B (K) -B (J)*F (K) )*SINF (UJ (L) ) 
DO  12  1=1,3 
GO  TO  (13, 14, 15), I 

13  J=2 
K=3 
GO  TO  12 

14  J=3 
K=l 
GO  TO  12 

15  J=l 
K=2 

12  H2(I)=B(J)*G(K)-G(J)*B(K) 
CALL  DRCTN  (B,  H2,  V(LX),  XLINE) 
DO  601  1=1,3 

601  A(I,LY)=XLINE(I) 
DO  113  1=1,3 

113  X(I,LY)=X(I,LX)+E(LY)*A(I,LY) 
TEM1=UJ(L)*57. 29577953 

PRINT  203,LX,X(1,LY),X(2,LY),X(3,LY),TEM1 
203  FORMAT  (15 , 3F15 . 10,F12. 7) 

114  DO  17  LJ=1,L 
TEM2=0.0 
DO  16  1=1,3 
TEM1=X  (I,L+2)  -X  (I  ,LJ) 

16  TEM2=TEM2+TEM1*TEM1 
S=SQRTF(TEM2) 
JI=LJ-1 
PRINT  206,JI,S 

206  FORMAT  (I8,F15.10) 
IF  (LJ-L+2)30,30,38 
30  IF  (S-R(L+2)-R(LJ))32,17,17 
32  PRINT  34 

34  FORMAT  (/6HREJECT//) 
IF  (SENSE  SWITCH  1)35,36 

35  PAUSE 

36  GO  TO  22 

17  CONTINUE 
38  IF  (SEN6E  SWITCH  1)18,19 

18  PAUSE 

19  IF  (EH(LX))119,43,119 

10 


Fortran  II  Program  (Cont. ) 


119  PSI=(V(LX)-3.141592654)/2.0 
CALL  DRCTN  (B,  H2,  PSI,  H3) 
THETA=HCH(LX)/2.0 
DO  42  KH=1,2 

CALL  DRCTN  (H3,  G,  THETA,  H4) 
DO  40  1=1,3 

40  XH(I)=X(I,LX)+EH(LX)*H4(I) 

PRINT   207,L,XH(1),XH(2),XH(3) 

207  FORMAT    (6H  HI553F15.10) 
IF    (THETA) 42, 43,  <+2 

42  THETA= -THETA 

43  IF  (L-N-l)20,21,20 

21  PRINT  205,ISOMR 
205  FORMAT  (14/) 

KK=N+3 

DO  39  J=1,KK 
LX=J-1 
39  PUNCH  208,ISOMR,LX,X(1,J),X(2,J),X(3,J),NSMPL 

208  FORMAT  (15 ,I3,3F15. 10,22X,I5) 
PRINT  209 

209  FORMAT (//) 
ISOMR=ISOMR+l 

20  L=L+LD 

IF  (LD)22,23,23 

23  IF  (L-N-l)24,24,25 

25  LD=-1 

GO  TO  20 

24  UJ(L)=UO(L) 
GO  TO  500 

22  IF  (L-2)100, 26,26 

26  UJ(L)=UJ(L)+DU(L) 

IF  (UJ(L) -6. 283185)27, 25, 25 

27  LD=1 

GO  TO  500 
100  PAUSE 

GO  TO  101 

END 
*1005 

SUBROUTINE  DRCTN (FROM,  TO , ANGLE ,XLINE) 

DIMENSION  FROM(3),TO(3),XLINE(3) 

DO  600  1=1,3 
600  XLINE(I)=FROM(I)*COSF(ANGLE)+TO(I)*SINF (ANGLE) 

RETURN 
END 
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Specifications  for  Fortran  Program 


(1)  The  program  computes  one  chain  of  atoms  having  individual  characteristics  with  two 
(or  less)  identical  side  atoms  symmetrically  attached  to  each  chain  atom  except  the 
first  and  last. 

(2)  Maximum  length  of  chain  =21  atoms.     Maximum  value  of  n  =  18. 

(3)  Number  of  side  atoms  per  chain  atom  is  governed  by  ^_HCpH  and  LHp.  If  both  are  non- 
zero, two  identical  side  atoms  will  be  placed  symmetrically  on  chain  atom  p.  If  /_HCp  H 
is  zero,  only  one  side  atom  will  be  attached,  as  for  a  carbonyl  oxygen  or  an  olefinic 
hydrogen.  If  LHp  is  zero,  no  side  atom  will  be  placed  on  the  chain  atom,  as  for  the 
chain  oxygen  atom  of  an  ether,    alcohol,  or  ester. 

(4)  Since  a  maximum  of  two  side  atoms  can  be  placed  on  a  chain  atom,  a  three-atom  set 
such  as  the  hydrogens  of  a  methyl,  at  one  end  of  the  chain,  requires  that  chain  atom 
numbered  zero  or  n+2  must  be  a  hydrogen. 

(5)  Input:  All  angles  in  degrees,  all  lengths  in  the  same  unit,  such  as  Angstroms,  in  order 
of  entry   shown  on    page  13. 

(6)  Output:  Printed  output  gives  coordinates  and  interatomic  distances  for  each  new  location 
of  each  chain  atom  and  coordinates  of  each  side  atom.  Each  location  is  printed  only 
once,  although  the  atom  at  this  location  may  be  involved  in  several  isomers.  Punched 
card  output  consists  of  a  complete  summary  of  the  coordinates  of  the  chain  atoms  for 
each  rotational  isomer,  for  use  in  further  computations  of  molecular  properties  such  as 
electric  moment,  interaction  energy,  moment  of  inertia,  etc.  Output  format  is  given 
below. 


Printout  format: 

(Chain  Atom  No.  ) 
(0) 
(1) 
(2) 


Output  Format  for  Fortran  Program 


(3  coordinates  of  numbered  atom  of  chain)  (Angle  J  of  chain) 

(Distance  between  numbered  atom  and  Atom  0) 
(Distance  between  numbered  atom  and  Atom  1) 
(Distance  between  numbered  atom  and  Atom  2) 


H    (Chain  Atom  No.  ) 
H    (Chain  Atom  No.  ) 


(3  coordinates  of  side  atom  attached  to  numbered  chain  atom) 
(3  coordinates  of  other  side  atom) 


(Isomer  No.  ) 
Punched  card  format: 
(Isomer  No.  )  (Chain  atom  No.  )    (3  coordinates  of  numbered  chain  atom)    (Sample  No.  ) 


12 


Card  Inputs  Required  by  Fortran  Program 

Using  the  Symbolism  of  the  Flow  Diagrams 

In  order  of  reading: 

n  Sample  No. 

xQ  y0  z0  Format   for   all  these   numbers     is 

F14.  10,    except  n,    which  is  12,    and 
sQ  tQ  uQ  Sample  No.  ,   which  is  14. 

x-^  yj  z^  Note  that  the  sequencing  places  only 

two  side  atoms  on  each  chain  atom. 

s^  tj  u^  Therefore,  in  order  to  accommodate 

a  methyl  group  at   each  end  of  the 

0.0  0.0  0.0  0.0  R0  chain,    the  chain  atoms  numbered  0 

and  n+2  must  be  hydrogens. 

0.0  0.0 

In  order  to  accommodate   a  single 

0.0  Vi  J\    o  AJi  R\  side  atom,    such  as  a  carbonyl  oxy- 

gen,   set  the     corresponding    HCH 

/HCiH  LHj  angle  =  0.0. 

L2  V2  J2    0  ^2  ^2  ^n  or<^er  t°  accommodate   a    chain 

atom  having  no   side   atoms,     such 

^/HCoH  LH9  as  an  ether  oxygen,    set  both/  HCH 

and  LH  =  0.  0. 

L3  V3  J3j  0  AJ3  R3 


/hc3h  lh3 


Ln 

vn 

/.HCnH 

LHn 

Ln+1 

Vn+1 

ZHCn+lH 

LHn+1 

Ln+2 

0.  0 

0.0 

0.0 

Jn,  0  AJn  Rn 


0.0  0.0  Rn+1 


0.0  0.0  Rn+2 
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Item 


Table  1 
Correlation  of  Systems  of  Symbols 

J.    Chem.    Phys.  Flow  Diagrams 


Fortran 


Serial  number  of  atom 

Maximum  value  of  p 

Maximum  value  of  L 

Serial  number  of  atom 
during  distance  loop 

Increment  or  decrement 

Coordinate  number 
subscript 

Direction  cosines  of 
bond  AB  pointing  from 
atom  p-1  to  atom  p. 

Direction  cosines  of 
bond  BC  pointing  from 
atom  p  to  atom  p+1 

Direction  cosines  of 
bond  CD  pointing  from 
atom  p+1  to  atom  p+2 

Direction  cosines 
of  line  BF 


Direction  cosines 
of  lines  BG  and  CG' 


Coordinates  of  atom  p+1 


Coordinates  of  atom  p+2 


0,    1.    2 


ao 
"1 

»2 

£o 
Pi 
£2 

?o 

Tl 

72 

*o 
*1 

*2 

uo 
Mi 

M2 

Co 

Ci 

C2 

Do 
Dl 

D2 


p  =  0,    1,    2.  .  .       L  =  p  (basic  comput.    unit) 
L  =  p+1  (cycling  and  input - 
output) 


n 

N  =n 

N+  1 

q=  0,    1,    2.. 

.       LJ  =  q  +  1 
JI    =  q 

Ap 

LD 

1=   1,    2,    3 

SP 
*P 

UP 

A(1,L) 
A(2,L) 
A(3,L) 

sp+l 
tp+1 

up+l 

A(l,    L+1)  = 
A(2,    L+1)  = 
A(3,   L+1)  = 

B(l) 
B(2) 
B(3) 

sp+2 
Up+2 

A(l,    L+2) 
A(2,    L+2) 
A(3,    L+2) 

*o 
*2 

F(l) 
F(2) 

F(3) 

"2 

G(l) 
G(2) 
G(3) 

xp+l 

yp+i 
xp+i 

X(l,    L+1) 
X(2,    L+1) 
X(3,    L+1) 

xp+2 

yp+2 

zp+2 

X(l,    L+2) 
X(2,    L+2) 
X(3,    L+2) 
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Table  1  (Continued) 
Correlation  of  Systems  of  Symbols 


Item 


Supplement  of 
bond  angle  at  p 

Supplement  of 
bond  angle  at  p+1 

Length  of  bond  between 
atom  p+1  and  atom  p+2 

Dihedral  angle 

Initial  value  of 
dihedral  angle 

Increment  for 
dihedral  angle 

Angle  between  bonds  to 
side  atoms  on  chain  atom  p 

Length  of  bond  between 
side  atom  and  chain  atom  p 


J.    Chem.    Phys. 


Flow  Diagrams 


V 


L 
J 


p+1 


Lp+2 


'po 


AJT 


'P 
/HCpH 

LH„ 


Table  2 
Coordinates  of  hydroxyl  group  for  selected  rotational  isomers 


Fortran 
V(L) 

V(L+1) 

E(L+2) 
UJ(L) 

UO(L) 

DU(L) 
HCH(L) 

EH(L) 


Hydrogen  atom 

Oxygen  atom 

X 

y 

z 

X 

y 

z 

-0.  760 

1.  754 

±1. 146 

-0.  205 

2.  537 

±1. 

168 

-0.  689 

1.  708 

±1.291 

-0.  166 

2.  512 

±1. 

248 

-0. 037 

1.  285 

±2. 069 

0.  522 

2.  065 

±2. 

068 

0.  085 

1.  206 

±2.  153 

0.  589 

2.  022 

±2. 

115 

0.976 

0.  627 

±2. 438 

1.  530 

1.  411 

±2. 

415 

1.  116 

0.  536 

±2.438 

1.  607 

1.  361 

±2. 

415 

Carbon  atoms 

of  the  double  bond 

-0. 6685 


0. 6885 
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